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SUMMARY 


The  work  done  during  the  first  quarter  on  "Research  on  Unsteady 
Stall  of  Axial  Flow  Compressors",  Contract  No,  DA-49-186-AM-1  3(X), 
for  the  Harry  Diamond  Laboratories,  Washington,  D.  C.  ,  by  Cornell 
Aeronautical  Laboratory,  Inc,  is  reported  here. 

Initial  effort  in  this  program  is  concentrated  on  continuation  of 
previous  research  at  CAL  concerning  (1)  an  experimental  study  of  the 
separation  of  a  laminar  boundary  layer  from  walls  moving  upstream  and 
downstream  in  a  steady  flow  using  a  shrouded  rotating  cylinder  in  a  wind 
tunnel,  and  (2)  an  analytical  study  of  compatibility  conditions  between 
circulation  disturbances  on  an  actuator  disc  in  a  mean  flow  with  swirl 
and  the  appearance  of  vorticity  downstream  of  the  disc. 

The  rotating  cylinder  apparatus  was  assembled  and  checked  out. 

Certain  modifications  to  the  equipment  were  made  to  obtain  the  performance 
required  for  accurate  boundary-layer  measurements.  Pressure  distri¬ 
butions  have  been  obtained  for  both  the  unshrouded  and  shrouded  rotating 
cylinder.  Preliminary  boundary-layer  data  were  obtained  with  small 
static  and  total  head  probes. 

The  analysis  completed  during  the  previous  research  for  a  mean  flow 
with  free -vortex-type  swirl  through  a  single-stage  compressor  was  reviewed, 
and  an  error  in  the  basic  differential  equations  discovered.  It  was  determined 
that  the  error  did  not  affect  the  results  obtained.  The  differential  equations 
for  a  similar  case,  but  with  wheel-type  swirl,  were  set  up.  As  yet,  a 
solution  has  not  been  obtained  for  this  case.  Results  were  obtained  for  a 
simplified  flow  which  is  an  approximation  to  the  flow  through  a  single-stage 
compressor  in  a  narrow  annulus.  In  this  case,  it  was  found  that,  within  the 
limitations  of  the  small  perturbation  theory  used,  the  only  possible  propa¬ 
gating  type  of  disturbance  in  circulation  on  the  actuator  disc  is  one  fixed  in 
and  rotating  with  the  blades. 
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II 


INTRODUCTION 


Rotating  stall  has  been  a  continuing  problem  for  the  designers  of 
axial -flow  turbojet  engines  since  the  inception  of  their  use  in  aircraft. 

For  the  compressor  designer,  the  major  problem  associated  with  rotating 
stall  has  been  one  of  blade  fatigue  failure  resulting  from  blade  resonance 
with  the  passage  of  stall  zones.  From  a  performance  point  of  view,  unsteady- 
stall  phenomena  may  place  severe  limitations  on  the  acceleration  or  decel¬ 
eration  of  a  jet  engine. 

Extensive  research  effort  has  been  devoted  to  this  problem  during  the 
past  two  decades.  Progress  was  made  in  that  a  degree  of  understanding  of 
the  phenomenon  was  attained;  however,  the  basic  problem  remained  unsolved 
when  lack  of  funds  reflecting  a  declining  interest  in  jet  engine  development 
brought  this  research  effort  nearly  to  a  halt  in  the  late  1950's.  In  particular, 
the  unsteady  boundary-layer  separation  phenomena  which,  presumably, 
underlie  rotating  stall  of  compressor  elements  continue  to  pose  problems 
of  fundamental  interest. 

Due  to  a  revival  of  interest  in  axial-flow  jet  engine  development,  this 
program  was  initiated  to  continue  a  methodical  approach  to  solutions  for  the 
problems  posed  by  unsteady  compressor  stall.  The  work  reported  here  is 
->  continuation  of  a  research  program  concerning  rotating  stall  in  axial-flow 
compressors,  carried  out  by  the  Cornell  Aeronautical  Laboratory  over  a 
three -year  period  ending  in  January  1959.  The  present  program  is  a  research 
effort  for  the  purpose  of  increasing  understanding  of  rotating  stall  and  other 
unsteady  surge  characteristics  of  axial  compressors.  Included  in  this  effort 
will  be  studies  of  the  aerodynamics  of  compressor  stages,  studies  of  unsteady 
sectional -lift  characteristics  of  single -blade  elements  near  stall,  and  studies 
of  possible  methods  of  alleviating  unsteady  surge  characteristics.  In  the 
current  year's  effort,  the  initial  emphasis  is  on  the  continuation  of  two  studies 
started  under  the  original  program:  (1)  a  theoretical  study  of  the  role  of 
entering  swirl  pattern  in  causing  rotating  stall  of  a  single  stage,  and  (2)  an 
experimental  study  of  the  separation  criteria  for  a  moving  surface  (shrouded 
rotating  cylinder). 
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Ill  WORK  ACCOMPLISHED 


A,  Rotating  Cylinder  Experiments 

This  work  is  a  continuation  of  the  research  reported  in  Reference  1. 

The  experiments  were  motivated  by  the  need  to  establish  criteria  for  the 
conditions  for  separation  in  a  quasi-steady  laminar  boundary  layer  in  the 
case  involving  a  moving  separation  point.  The  theoretical  work  leading  up 
to  the  experimental  program  was  reported  in  Reference  2. 

During  the  previous  work,  it  was  found  in  analyzing  the  quasi -steady 
boundary  layer  near  separation  that  the  problem  involving  a  moving  separation 
point  could  be  reduced  to  a  steady-flow  problem  in  which  the  wall  is  in  motion. 
This  problem  was  treated  separately  for  the  cases  of  the  wall  in  downstream 
motion  and  in  upstream  motion.  In  the  former  case,  Sketch  A,  it  was  found 
that  vanishing  shear  occurred  at  some  distance,  proportional  to  wall  velocity, 
upstream  of  separation,  and  separation  was  characterized  by  a  separated 
region  embedded  in  the  tlow. 


SEPARATED 

REGION 


WALL 

MOVEMENT 


SKETCH  A 

WALL  MOVING  DOWNSTREAM 
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In  this  instance,  separation  would  not  occur  at  the  wall,  but  at  a  finite  height 
above  the  wall,  and  there  would  be  a  boundary -layer  flow  beneath  the  wake 
region.  The  analysis  provides  quantitative  estimates  of  the  displacement  of 
the  separation  point  due  to  wall  velocity  and  of  the  position  of  the  separation 
point  above  the  wall.  These  results  depend  somewhat  on  the  assumed  velocity 
profile  at  separation,  and  are  valid  only  for  wall  velocities  which  are  small 
in  comparison  with  the  free -stream  velocity.  Accordingly,  the  experimental 
program  was  set  up  to  determine  the  validity  of  the  assumed  profile,  and  to 
measure  the  effects  of  small  and  large  wall  velocities  on  the  position  of 
separation  as  a  check  on  the  theory. 

In  the  case  of  the  upstream  moving  wall,  Sketch  B,  it  was  hypothesized 
that  a  probable  profile  at  separation  was  one  that  flattened  against  the  axis 
over  a  relatively  large  height  with  a  singularity  at  the  foot;  that  is,  one  in 
which  the  velocity  abruptly  changed  from  zero  near  the  wall  to  the  wall 
velocity  at  the  wall. 


WALL 

MOVEMENT 


SKETCH  B 

WALL  MOVING  UPSTREAM 

In  this  model,  the  boundary  layer  then  contains  a  sublayer  of  circulating  flow 
extending  from  the  forward  stagnation  point  to  the  separation  point  and, 
rigorously,  there  would  be  no  wake  flow  into  the  boundary  layer.  The  theory 
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indicates  that,  in  this  case,  the  separation  point  would  be  displaced  upstream 
a  distance  proportional  to  the  wall  velocity,  but  provides  no  quantitative 
estimate  of  the  displacement.  It  is  noted  that  this  result  is  somewhat  specu¬ 
lative  in  that  it  is  largely  dependent  on  the  assumed  profile  at  separation.  In 
order  to  clarify  the  upstream  moving  case,  the  experiments  were  designed  to 
investigate  the  velocity  profile  near  separation  to  assess  the  validity  of  the 
assumed  singular  profile.  It  was  anticipated,  of  course,  that  a  singular 
profile  in  the  mathematical  sense  would  not  be  observed  and  that  there  would 
be  some  flow  from  the  wake  into  the  boundary  layer.  If  the  height  of  this 
reverse  flow  were  found  to  be  small  in  comparison  with  the  height  over  which 
the  profile  was  flattened,  then  the  assumed  profile  would  be  substantially 
verified.  If  this  were  not  the  case,  the  experimental  profiles  would  serve 
as  a  basis  for  improving  the  analysis.  Also,  since  the  analysis  provided  no 
quantitative  estimates  of  the  displacement  of  the  separation  point  due  to  wall 
velocity,  the  experiments  were  to  provide  this  information. 

Since  the  analysis  demonstrated  the  equivalence  of  the  quasi -steady 
problem  and  the  problem  of  a  moving  wall  at  separation,  the  model  chosen 
to  experimentally  investigate  unsteady  separation  was  a  rotating  cylinder  in 
steady  flow.  This  model  simulates  the  pertinent  features  of  the  unsteady 
problem  in  that  the  wall  and  separation  point  are  in  relative  motion,  and 
both  upstream  and  downstream  moving  walls  are  present.  Of  course,  the 
chief  feature  of  this  experiment  is  that  it  is  a  steady-flow  experiment,  and 
more  precise  measurements  are  feasible.  The  method  suffers  somewhat  in 
that  a  relatively  steep  pressure  gradient  is  present  at  separation,  and  the 
effects  of  a  moving  wall  on  delaying  or  promoting  separation  might  be  difficult 
to  evaluate  if  they  were  small  in  comparison  with  the  effect  of  the  strong 
pressure  gradient.  Also,  because  of  the  Magnus  effect,  the  pressure  distri¬ 
bution  depends  on  the  wall  velocity,  and  observed  shifts  in  separation  might 
be  th^  result  of  a  change  in  pressure  distribution.  In  order  to  circumvent 
these  difficulties,  provisions  were  made  to  shroud  the  model  so  that  a  more 
moderate  airfoil-type  pressure  gradient  could  be  imposed  on  the  cylinder. 

This  would  not  necessarily  remove  the  Magnus  effect,  but  it  was  anticipated 
that,  by  measuring  the  pressure  distribution  on  the  model  as  a  function  of 
rotation,  it  would  be  possible  to  change  the  shroud  shape  and  maintain  a 
given  pressure  distribution. 
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Experimental  Apparatus 


The  basic  apparatus  is  the  equipment  that  was  used  for  the  previous 
investigation  and  consists  of  a  4 -in.  diameter  cylinder  mounted  between 
plastic  endplates,  spanning  the  full  width  of  the  subsonic  test  section  of  the 
CAL  One-Foot,  High-Speed  Wind  Tunnel.  The  cylinder  is  mounted  on  bearings 
so  as  to  rotate  about  its  centerline. 

The  cylinder  is  instrumented  with  a  piezoelectric  pressure  transducer 
vented  through  a  1/8-in.  diameter  hole  in  the  surface  of  the  cylinder.  The 
usable  frequency  range  of  the  pressure  transducer  and  vent  are  0  -  40,  000  cps; 
however,  filters  are  incorporated  in  the  recording  system  to  limit  response 
to  the  range  between  one  cps  and  1000  cps  in  order  to  eliminate  steady-state 
and  slowly  changing  centrifugal  effects,  and  spurious  electrical  and  acoustic 
noise.  The  output  of  the  pressure  transducer  is  displayed  on  a  Tektronix 
Type  502  Oscilloscope  and  is  recorded  photographically.  Rotation  speeds  of 
the  cylinder  are  measured  by  means  of  a  miniature  electromagnetic  pickup 
in  conjunction  with  a  60 -tooth  gear  attached  to  the  cylinder  shaft.  The  output 
of  the  pickup  is  displayed  on  the  oscilloscope  as  well  as  on  a  Berkeley  Counter 
and  Timer. 

A  substantial  portion  of  the  work  of  this  first  quarter  has  been  devoted 
to  reassembly  and  checkout  of  the  apparatus.  Preliminary  tests  of  the  rotating 
cylinder  showed  that  the  performance  of  the  Onsrud  Air  Turbine  cylinder  drive 
motor  would  be  unsatisfactory  for  boundary  layer  survey  work,  where  constant 
rpm  must  be  held  closely  for  extended  lengths  of  time.  The  rotational  speeds 
varied  continuously  over  a  wide  range.  Furthermore,  fine  control  of  rpm  was 
extremely  difficult.  Therefore,  the  air  turbine  was  replaced  by  a  Heller 
Electrical  Motor  and  Controller  Unit,  Model  2T60.  Subsequent  tests  have 
shown  that  drift  of  this  unit  is  very  small,  and  fine  control  is  possible. 
Rotational  speed  variations  can  be  kept  within  +  1%  without  difficulty. 

Also  during  the  preliminary  shakedown  tests,  it  was  found  that  the 
cylinder  driveshaft  was  bent  slightly.  This  caused  excessive  vibration  and 
runout,  or  eccentricity,  of  the  cylinder,  A  new  driveshaft  was  fabricated,  and 
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one  of  the  bearings  replaced.  Subsequent  tests  showed  that  the  vibration  had 
been  reduced  to  a  satisfactory  level.  Also,  runout  of  the  cylinder  surface  is 
now  a  maximum  of  0.  0003  in.  ,  which  is  quite  satisfactory. 

Shrouds  were  available  from  the  previous  program,  as  noted  above. 

These  shrouds  are  mounted  above  and  below  the  cylinder  and  were  designed 
so  as  to  impose  a  more  gradual,  airfoil-type  pressure  gr  ent  than  exists 
on  the  unshrouded  cylinder. 

It  is  intended  that  the  bulk  of  the  boundary  layer  velocity  survey  work 
will  be  done  with  hot-wire  anemometer  equipment.  This  type  of  velocity  probe 
is  required,  rather  than  the  conventional  total  pressure  probe,  because  the 
essential  boundary  layer  velocity  data  must  be  obtained  in  the  approximately 
0.  01  in.  layer  immediately  adjacent  to  the  cylinder.  However,  small  diameter 
total  head  and  static  probes  have  been  adapted  to  the  traversing  mechanism* 
and  are  also  being  used  in  the  boundary  layer  investigations.  The  total  head 
probe  is  particularly  useful  in  defining  the  outer  edge  of  the  boundary  layer, 

The  flow  velocity  outside  the  boundary  layer  may  vary  continuously  with 
distance  away  from  the  cylinder  surface,  but  the  total  pressure  approaches 
a  constant  value  as  the  outer  edge  of  the  boundary  layer  is  reached. 

Details  of  the  CAL  One -Foot,  High-Speed  Wind  Tunnel  are  presented 
in  References  1  and  3.  The  configuration  of  the  wind  tunnel  being  used  in  the 
current  tests  are  as  described  in  Reference  1;  stream  turbulence  level,  and 
test  section  calibrations  indicated  in  Reference  1  should  be  applicable  to  the 
present  work. 

Complete  details  of  the  experimental  apparatus  are  included  in  Reference 
1.  Those  minor  modifications  which  have  been  incorporated  during  the  present 
program  have  been  noted  above. 


*  The  traversing  mechanism  is  a  jeweler's  lathe  slide -rest,  or  compound, 
mounted  on  top  of  the  wind  tunnel  test  section.  Slots  cut  in  the  top  of  the 
test  section  and  in  the  top  shroud  permit  the  probes  to  be  used  in  the 
proximity  of  the  cylinder  surface. 
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Experimental  Data  Obtained  to  Date 


Figure  1  illustrates  the  notation  used  for  the  rotating  cylinder  experi¬ 
ments.  With  the  free-stream  velocity  as  shown,  cylinder  angular  rotation 
velocity  is  positive  counterclockwise;  the  angle  on  the  cylinder  is  measured 
from  the  horizontal  reference  (radius  parallel  to  free-stream  direction), 
positive  in  a  counterclockwise  direction;  angles  between  0°  and +180°  are 
on  the  upper  surface  of  the  cylinder,  and  between  0°  and  -180°  on  the  lower 
surface.  In  the  boundary  layer,  y  is  the  distance  away  from  the  surface  of 
the  cylinder  along  the  local  normal  to  the  surface,  V  is  the  peripheral  speed 
of  the  surface  of  the  cylinder,  and  Vb  is  the  local  velocity  in  the  boundary 
layer  parallel  to  the  tangent  plane  to  the  cylinder  surface.  U  is  the  free- 
stream  velocity. 

Pressure  distributions  have  been  obtained  for  the  unshrouded  rotating 
cylinder  at  various  rpm's.  These  data  are  essentially  in  agreement  with 
corresponding  data  presented  in  Reference  1.  The  unsteady  pressure 
fluctuations  noted  in  Reference  1  for  the  unshrouded  rotating  cylinder, 
attributed  to  vortex  shedding,  were  also  obtained. 

Pressure  distributions  for  the  shrouded  rotating  cylinder  have  been 
obtained  over  a  range  of  positive  and  negative  rpm.  This  type  of  data  had  not 
been  obtained  during  the  previous  work.  Figure  2  shows  a  typical  oscillograph 
record  of  traces  from  the  output  of  the  cylinder  pressure  transducer  for 
-O.  =  -1284  rpm,  V  =  58.  3  fps;  this  corresponds  to  a  V  /  ZT  of  0.  385. 

Also  shown  in  Figure  2  are  the  corresponding  wind -off  calibration  traces.  In 
both  cases,  the  traces  shown  include  at  least  five  consecutive  double  cycles 
superimposed  one  on  top  of  the  other;  hence,  the  results  shown  in  Figure  2 
are  indicative  of  repeatability  of  the  pressure  data  and,  also,  steadiness  of 
the  flow  over  the  shrouded  cylinder.  Also  apparent  in  the  data  of  Figure  2  is 
a  certain  amount  of  contact  noise,  evidenced  by  the  relatively  high  frequency 
uneveness  in  the  traces.  Note,  in  particular,  the  large  "pip"  in  the  middle 
of  both  sets  of  traces  in  Figure  2.  The  appearance  of  such  contact  noise  in 
the  traces  has  some  consistency,  but  not  enough  to  be  accurately  calibrated 
out.  Therefore,  when  this  data  is  reduced,  a  smooth  faired  curve  is  drawn 
through  what  appears  to  be  the  average  of  the  traces  recorded. 
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Figure  3  shows  the  static  pressure  distribution  in  terms  of  the  pressure 
coefficient  for  the  same  conditions  as  the  Figure  2  data.  Here,  &ps  *s  the 
difference  between  the  local  and  the  free -stream  static  pressures,  and  q  is 
the  free-stream  dynamic  pressure.  The  same  technique  is  being  used  to  reduce 
the  oscillograph  data  as  was  used  in  Reference  1.  The  plotted  data  in  Figure  3 
was  obtained  from  averages  of  at  least  five  traces  for  two  successive  cycles 
in  each  trace.  The  differences  between  the  two  cycles  of  data  are  believed 
representative  of  scatter  resulting  primarily  from  the  contact  noise  and  the 
data  reduction  method  which  must  be  used;  these  data  show  the  largest  scatter 
obtained  from  the  shrouded  cylinder  data  thus  far  reduced.  Figure  4  shows 
similar  data  for  TJ  =  58.  3  fps,  _TL  =  +712  rpm,  where  much  better 
agreement  between  the  time -averaged  data  for  two  cycles  was  obtained.  Also 
shown  on  Figures  3  and  4  are  static  pressure  coefficient  data  obtained  with  the 
static  probe  at  Q  -  +  75°.  Agreement  is  very  good  for  fl  =  +712  rpm, 
and  within  the  scatter  for  iT  =  -1284  rpm. 

The  shrouds  appear  to  have  produced  the  more  gradual  pressure 
gradient  desired.  The  repeatability  evidenced  in  the  photographic  data 
appears  to  confirm  that  gross  unsteadiness  in  the  flow  due  to  vortex  shedding 
has  been  surpressed. 

There  does  appear  to  be  some  variation  in  the  pressure  distribution 
with  rpm.  Whether  this  change  in  pressure  distribution  itself  can  account 
for  the  shift  in  separation  point  with  change  in  rotation  speed,  for  which  data 
was  presented  in  Reference  1,  will  require  further  investigation.  It  should  be 
remembered  that  the  intent  of  this  experimental  program  was  not  only  to 
investigate  the  boundary  layer  velocity  profile  at  separation  for  upstream  and 
downstream  motion  of  a  separation  point,  but  also  to  check  the  results  of  the 
theory  of  Reference  2  as  regards  the  relative  motion  of  the  separation  point 
itself  with  wall  velocity  variation  for  walls  in  motion  downstream,  and  to 
obtain  similar  data  for  the  upstream  case  where  the  analysis  provided  no 
quantitative  estimates. 
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Figure  t  Apjfvs  6  FOR  SHROUDED  ROTATING  CYLINDER 
SI  =  712  RPM  (Re)p  =  1.2  x  I05 


Velocity  profiles  have  also  been  obtained  for  the  outer  portion  of  the 
boundary  layer  at  &  =  90°  for  a  range  of  cylinder  rotational  speeds  between 

H  =  0  and  /l  =  -1750  rpm.  This  data  is  not  presented  here  because  it  is 
incomplete;  the  probe  cannot  accurately  measure  the  flow  in  the  critical 
region  of  the  boundary  layer  near  the  cylinder.  For  this  region,  hot-wire 
equipment  must  be  used.  However  the  total  head  probe  data  will  be  useful 
as  a  check  on  the  hot-wire  results  and  will  be  particularly  useful  in  defining 
the  outer  edge  of  the  boundary  layer,  where  the  total  pressure  approaches 
the  constant  value  of  the  outer  flow,  regardless  of  the  variation  in  velocity 
in  the  outer  flow. 


B.  Theoretical  Analyses 

The  initial  theoretical  effort  during  this  program  is  concerned  with  a 
continuation  of  the  work  reported  in  Reference  4.  This  work  was  a  study  of 
possible  disturbance  patterns  which  may  exist  in  the  flow  field  produced  by 
a  single  compressor  stage  — with  a  view  to  the  description  of  possible  three- 
dimensional  or  partial -span  modes. 

It  should  be  made  clear  that  this  theoretical  study  is  not  directed  toward 
attempts  to  predict  specific  characteristics  of  rotating  stall  phenomena,  such 
as  speed  of  propagation,  number  of  stall  zones,  or  the  prediction  of  when 
rotating  stall  will  occur  in  terms  of  compressor  characteristics.  Rather,  the 
investigation  is  concerned  with  compatibility  conditions  in  the  flow  which 
determine  whether  rotating  stall  is  a  permissible  disturbance  of  the  particular 
type  of  flow  considered. 

It  was  found  in  Reference  4  that,  in  a  two-dimensional  cascade,  stall 
must  propagate  normal  to  the  blade  elements  and  cannot  propagate  in  an 
oblique  direction.  For  the  single-stage  axisymmetric  case  having  a  mean 
swirl  of  the  free -vortex  type,  it  was  found  that  no  partial -span  disturbance 
of  the  rotating  stall  type  was  possible  and,  further,  that  only  disturbances 
propagating  with  the  blade  rotation  velocity  were  possible.  These  are  potential 
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flow  disturbances  which  are  contingent  on  a  combination  of  blade  element 
negative  lift-curve  slope,  and  steady  hysteresis  of  the  lift  curve. 

Initial  effort  during  the  first  quarter  was  devoted  to  an  attempt  to 
obtain  solutions  to  the  basic  differential  equations  treated  in  Reference  4 
in  which  the  assumed  mean  flow  included  wheel-type  swirl,  rather  than  free- 
vortex-type  swirl  as  treated  in  Reference  4.  During  the  early  stage  of  this 
work,  the  differential  equations  reported  in  Reference  4  were  checked,  and 
it  was  found  that  there  was  an  error  in  one  of  the  terms.  It  was  felt  that, 
because  of  this,  there  could  be  a  change  in  the  results  of  the  free -vortex-type 
swirl  analysis;  therefore,  the  analysis  was  repeated  based  on  the  corrected 
equations 


Corrected  Free -Vortex  Swirl  Analysis 

The  notation  and  coordinate  systems  used  here  are  the  same  as  those 
of  Reference  4.  Cylindrical  coordinates  (  %>  ,  r  ,  Q  )  are  used  with  the 
X-  -axis  centered  in  the  compressor  disc.  The  mean  axial  flow  velocity 
component  is  U0  ,  and  there  is,  in  addition,  a  mean  tangential  flow  velocity 
U0S(r)  .  The  compressor  disc  is  represented  by  an  axisymmetric  actuator 
located  at  X>  =  0. 

The  general  equations  for  the  convection  of  vorticity  (Eq.  27  of  Reference 
4)  are: 

-J~  +(?'V)&  «  (OJ  -V)%  (1) 

where  W  is  the  vorticity  vector  and  ^  is  the  velocity  vector.  We  write 

fS  U0  ■*  U  ,  IT ,  U0S(r)  +  ur  (2) 

UJ  a  ^-(rS)r  +  S,q,C  (3) 

where  (  U  ,  If  ,  UT )  and  (  f  ,  f)  ,  £  )  are  perturbation  components,  and 
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(4) 


£  -j(ru/)r  -  7  ir6 

1 

7*7  “a  ~ 

*  “  -  OLr 

from  the  relation 

CL  *  V  *  ^ 

Equation  (1)  becomes,  after  linearization, 


i 

(A  +  u  <9  ,  UoS  9  \ 
\dt  U°3X  +  r*  39) 

7 

■  < 

^-(rS)rU', 

f. 

~r(^)f  v-  —-CrS^  7  -  2-^  IJ 

which  corresponds  to  Eq.  (33)  of  Reference  4;  the  underlined  term  on  the 
right-hand  side  in  the  equation  for  is  a  correction.  For  vortex  swirl, 


S* 


K  a  constant,  Eqs,  (5)  become 


(JL 

\dt 


u0 


dX. 


U0X  d_\ 

r 2  del 


n 


(6) 


which  is  the  corrected  form  of  Eqs.  (35)  of  Reference  4.  The  solutions  for  jjj 
and  ti  are 

£  -  f(\,r) 

.  (7) 

7  -  y  (a  .  ^ 

as  before  (Eqs.  (36),  (37)  and  (40)  of  Reference  4)  where 


A -e-(~*-$[)x  +  +  u.e  (8) 

However,  we  must  replace  Eq.  (38)  of  Reference  4  for  the  £  component  of 
perturbed  vorticity  by  the  solution  of 


UoK  ± 

r‘  39 


(9) 
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subject  to 


Ct-T*+ 


(10) 


The  "+"  subscript  refers  to  conditions  just  downstream  of  the  compressor 
disc,  i.  e.  ,  at  X  =  0+.  A  solution  of  Eq.  (9)  subject  to  Eq.  (10)  is 

(ID 

which  replaces  Eq,  (38)  of  Reference  4. 

Equation  (43)  of  Reference  4  now  becomes 

-0  (12) 

which,  at  the  actuator  disc  (X  =  0)  is 

+  -0  (13) 

as  before.  The  remainder  of  the  analysis  proceeds  as  before;  the  essential 
change  in  the  results  reported  in  Reference  4  for  free-vortex  swirl  is  the 
added  second  term  on  the  right-hand  side  of  Eq.  (11),  which  disappears  at  the 
actuator  disc  and,  hence,  does  not  affect  the  succeeding  analysis  which  is 
based  on  conditions  at  the  disc. 


Wheel-Type  Swirl  Analysis 
For  wheel-type  swirl, 

Sfr)  -  <r  (14) 

and  Eqs.  (5)  become 


(  d  d 

U  +yVx 

+  KU°<?9  )' 

’I' 

.  - 1<. 

’vl/  ’ 

^  1 

-C. 
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Equations  (15)  plus  Eq.  (4)  relating  (  £  ,  q  ,  £  )  to  (  u  ,  if  ,  w  ) 
constitute  the  six  simultaneous  differential  equations  relating  the  three 
perturbation  vorticity  components  and  the  three  perturbation  velocity 
components  for  wheel -type  swirl. 


If  the  independent  variables  are  transformed  as 


* 


r  « 


uet 

R0 


and  if  . 

u.  -  Do  it.' 

V  ■  U0  ir‘ 

to  -  U„uj‘  , 

on  substituting  Eq.  (4)  in  Eqs.  (15),  we  obtain 


_  dif'l 


DTidr,(ru/^  ae  J 


2Kr 


p  r  3 

DZ 


/  a  it! 

dr! 

,  dir' 


D  du!  dir'  ~\  /  a  Hr' 

DZ  L dr'  ~  &K'  J  ‘  diu' 


where 


_D_ 

DZ 


«  A  +  jL  +  <A. 
~~az  dx'  ^dQ 


(16) 


It  is  possible  to  derive  from  Eq.  ( 1 6)  partial  differential  equations,  each 
in  terms  of  U>  ,  If'  and  lO'  only,  but  of  higher  order.  That  for  U.'  is 


D  [  d  [  ,  deu'  r'(r'-f  1)  33 u  1) 

Dr}dr'[r  dr'ae  r'-t  dx'*de\] 


(rV 

r'-l 


iluLI 

Bx,l\ 


=  0 


(17) 
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and  that  for  ir '  is 


D  f  2r'  d^(r'is')  l-r'  d+(n'v')  1  .jJVV) 

DT  \r'+1  dvi'dr'  r'(r'+1)  dx' dr'dQ1-]  *  di.'*d&  * 


(18) 


None  of  the  sets  of  equations,  (4)  and  (15),  ( 1 6),  or  (17)  and  (18),  have  as  yet 
been  solved,  nor  does  it  appear  that  solutions  will  be  obtainable  without 
simplifications. 


Analysis  for  Mean  Flow  with  Swirl  Through  a  Narrow  Annulus 

If  it  is  assumed  in  Eq.  (5)  that  5  m  K  (a  constant),  V -  O  and  that  u 
and  us  do  not  vary  with  f  ,  then  the  flow  could  be  considered  as  approxi¬ 
mately  that  through  a  very  narrow  annulus.  If  this  annulus  is  centered  at 
R0  ,  Eqs.  (5)  can  be  written 


(-2-  +  U  Z  +  UdLZ) 

\at  u‘Sx  R,J9l 


21 

du  i 

dx 

7 

KU0 
Zo  . 

0 

dus  n 

.  C. 

-n) 

As  from  Eq.  (4) 


d  us  1  d\s 

av'  '  ae  = 


1  du 

"  ~  R.0  d& 


du. r 

<?Z 


4  ~  ax  "  dr 


we  have  from  Eq.  (19), 
dU 

* 

or 


3  us  _ 

ax> 


7  * 


(19) 


(20) 


(21) 


U,  - 


(0,-b) 


(22) 
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and  from  Eqs.  (20)  and  (21) 

duJ 


(23) 


_  _  1  du 

dv  “  2R0  99 


The  continuity  equation  is 


+  ±  ii"  x  0 
9X  R.09Q 


from  which,  as 


(24) 


or 


its  -  us  (x,i) 


(25) 


From  Eq.  (19) 


+  U0 


9g  Do  K  9q 

dx  +  K0  Te 


o 


so  that 

7  *  0 (1)  » 

where 


A  =  9 


-M) 


X  4-  a  0o  i 


(26) 


As  in  Reference  4,  "f  is  defined  as  the  circulation  per  unit  angle  on  the 
actuator  disc.  For  a  specific  blade  element  (Eq.  45  of  Reference  4) 

-  "•*•?*  <27> 

where  kUo  is  the  rotational  velocity  of  the  blade  elements. 


Let 

7(9. t)  -  7'(\*.  >  c r) 

where 

X+  -  0+*Uet 
i 7"  •  &  -  kVo  4 
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Under  this  transformation,  which  is  permissible  except  for  -a-  =  -  k  ,  Eq.  (27) 
becomes 


dr' 

2\+ 


4,  +  k 


9^t) 


(28) 


which  can  be  integrated  to  give 

y'  ~  ~~"k  J  g(v)  dz -h  m(cr) 

We  have  from  Eq.  (25)  that  ur  *  11/ (iC> 4)  and  from  Eq.  (21)  that  q 
or  that 

f)  =?(*)  = 


dur(.%,-i) 

3% 


(29) 


2  VO- 


(30) 


Here, 


hence, 


e 


can  at  most  be  a  constant, 


say,  g0 


du. 

36 


From  Eqs.  (23)  and  (30) 


u  = 

However,  as  u(6,t)mu(&  +  2Tr>4’)  ,  q0  must  be  zero. 

Thus, 

-f  -  m(<r)  =  m  (&  - kOet ) 

and  the  only  circulation  disturbance  of  the  rotating  stall  type  possible  is  one 
fixed  in  the  blades. 
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IV  WORK  FOR  THE  NEXT  QUARTER 


Experimental  Program 


Work  will  continue  with  the  rotating  cylinder  to  obtain  boundary  layer 
velocity  profiles  at  separation  for  the  wall  moving  both  downstream  and 
upstream.  The  development  of  the  boundary  layer  velocity  profile  as  the 
separation  point  is  approached  will  also  be  investigated.  Hot-wire 
anemometer  equipment  will  be  available  for  this  work. 


Theoretical  Program 

Work  will  continue  on  the  compatibility  conditions  relating  disturbances 
in  circulation  in  the  actuator  disc  to  the  appearance  of  vorticity  downstream. 
Initially,  the  narrow  annulus  configuration  will  be  emphasized;  attempts  will 
be  made  to  relax  the  assumptions  made  thus  far  (V  =  0,  no  variation  of  U. 
and  "W  with  radius)  which  are  undoubtedly  oversimplified. 
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